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ABSTRACT 

Accurate distances to pulsars can be used for a variety of studies of the Galaxy 
and its electron content. However, most distance measures to pulsars have been 
derived from the absorption (or lack thereof) of pulsar emission by Galactic 
H I gas, which typically implies that only upper or lower limits on the pulsar 
distance are available. We present a critical analysis of all measured H I distance 
limits to pulsars and other neutron stars, and translate these limits into actual 
distance estimates through a likelihood analysis that simultaneously corrects for 
statistical biases. We also apply this analysis to parallax measurements of pulsars 
in order to obtain accurate distance estimates and find that the parallax and H I 
distance measurements are biased in different ways, because of differences in the 
sampled populations. Parallax measurements typically underestimate a pulsar's 
distance because of the limited distance to which this technique works and the 
consequential strong effect of the Galactic pulsar distribution (i.e. the original 
Lutz-Kelker bias), in H I distance limits, however, the luminosity bias dominates 
the Lutz-Kelker effect, leading to overestimated distances because the bright 
pulsars on which this technique is applicable are more likely to be nearby given 
their brightness. 

Subject headings: astrometry, pulsars 
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Introduction 



The rotation of pulsars, which causes their continuous emission to be observed as 
highly regular pulses, makes these objects highly useful probes of any dispersive phenomena 
in interstellar space. Combined with an accurate and precise distance, pulsar emission 
(specifically its dispersion and Faraday rotation) provides crucial information for modelling 
of the Galactic electron distribution and magnetic field. 

Parallax measurem ents are non-trivial undertakings and only very few significant 



parallax measurements (jGwinn et al. 



1986 



Bailesetal 



19901 ) were made within the first 



two decades after pulsars were discovered. Another method to determine a pulsar's distance 
is based on Galactic H I spectra in the direction to the pulsar. This method (known as the 
kinematic or H I method) compares the H I spectrum on-pulse (when the pulsar emission is 
seen) and off-pulse (when the pulsar emission beam is turned away). Any observed pulsar 
absorption must originate in gas lying closer than the pulsar; while gas located farther 
than the pulsar will not exhibit absorption. The velocities of these respective H I regions 
are subsequently derived from the spectrum and translated to distances with help of a 
Galactic rotation model. The distance of the furthest H I gas that appears in absorption 
then provides a lower limit D\ ow on the pulsar distance, while the distance of the nearest gas 
that only appears in emission, is interpreted as an upper limit D np on the pulsar distance. 



Roughly two decades after the discovery of pulsars, iFrail fc Weisbergi (jl990j, henceforth 
FW90) collated all published pulsar distances, which at the time consisted of 50 H I 
distances, three parallax measurements and 20 distances by association. Given the 
importance of H I distances, they critically investigated the various measurements and 
defined a set of criteria that has been used in almost all subsequent publications. 

Progress in both interferometric hardware (at the Long Baseline Array in the South 
and the Very Long Baseline Array in the North) and in the sensitivity of pulsar timing, 



-4- 



subsequently allowed an exponential increase in the nu mber of measur e d pul sar parallaxes 



so that currently 57 parallaxes are measured. This led 



Verbiest et al. 



(120 lOl henceforth 



VLM10) to co ll ate t hose distances and invest igate the statistical bias predicted by 



VLMlOl was based on a Bayesian analysis 



Lutz fc Kelkeri (119731 ). The work presented by 
that took into account both the Galactic distribution of pulsars (which is the actual bias 
first discussed by Lutz and Kelker in 1973) and the intrinsic pulsar luminosity distribution; 
but they only considered parallax measurements. 



FW9d we list all 80 published 



In this paper, we present an update of the work done by 
distances to pulsars and other neutron stars, based on H I measurements or associations 
wit h obje cts having H I distances, and evalu ate them based largely on the criteria laid out 



by 



FW90 



. We then improve the analysis of 



VLMlOl by deriving fully analyti c solution s 



that replace the need for (approximate) Monte-Carlo simulations. Also, the 



VLM10 



analysis is expanded to incorporate information provided by H I distance lim its; and to 



VLMld 



provide bias-corrected distances in addition to parallaxes. As in the case of 
the present paper bases its bias-correction method on empirical models for the Galactic 
pulsar distribution and luminosity function. These models do add an unquantified level of 
uncertainty to the analysis, but can easily be updated as our knowledge about the pulsar 
population grows through pulsar surveys. The evaluation of H I distance limits is presented 
in Section [2j the likelihood analysis to correct for the biases is derived in Section |3j 
Bias-corrected parallaxes and distances are given in tables [1] and [2] and a summarising 
discussion is found in Section HJ 
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H I Kinematic Distances 



2.1. Source Selection and H I Kinematic Distance Limit Determination 



FW90I established standard techniques for the extraction of reliable pulsar kinematic 



distance limits. Specifically, they defined the bound of > 35 K on the brightness 
temperature of H I emission used for deriving upper distance limitJ3, pointing out that 
weaker emission would not be expected to result in significant absorption. Secondly, 
they re-evaluated distance limits based on old Galactic models and rotation curves, 
defaulting to the IAU values for the distance of the Solar System to the centre of the 
Galaxy [Rq = 8.5 kpc| and the Galactic rota tion velocity in the Solar System nei ghborhood 



[ 8n = 220 km s" 1 ] flKerr fc Lvnden-Bell 



19861 ). and using the flat rotation curve of 



FichetaL 



( 119891 ). In converting velocities to distances, furthermore, the y assumed an uncertainty 
of 7 km s -1 because of known random motions of that order ( IDickey fc Lockmanlll990l ). 
Finally, in the Perseus arm, with its well-known spiral shoc k, they either used i ndependent 



Joncas et al. 



f|l989h . which 



distance tracers or applied the approximation proposed by 
states the global rotation curve can be applied (near G\ = 130°) provided the measured H I 
velocities are decreased by a factor of 1.6. Most investigations since then have used these 
same criteria and so does the present paper, with a few exceptions as listed below. 

In the current work, we present a uniformly-determined sample of neutron star H I 
kinematic upper and lo wer dis tance estimates by finding all such efforts in the literature, 



and then applying the 



FW90 



criteria to any published data that have not previously 
been analyzed with that procedure. If the cited authors made a good case for a non-flat 
rotation curve (e.g, in the direction of the Galactic bar, the 3 kpc arm, or the Perseus arm 



x They also note this bound can be relaxed depending on the sensitivity of the observation, 
provided the optical depth for the emission is 0.3 or higher. 



-6 - 



shock), we maintain their curve in our analysis. If, however, the original authors used a flat 
rotation curve but non-IAU Galactic constants, we reanalyze t he kinemati c dista nce limits, 
using the flat rotation curve and IAU constants. We note that 



Reid et al. 



((20091) obtain a 



Galactic rotation velocity that is larger than the IAU value at a signficance of 95%. More 
recent measurements by the same authors have increased the signific ance of this offset t o 



close to 99% (Reid, personal communication). Based on Eq. 2.21 of iDamour fc Taylor 



(119911 ). we find that this could imply an overestimate of our H I distance limits by up to 
~ 20% (though generally much less), depending on Galactic longitude and the measured 
H I velocity. 

In Table [TJ we list the H I kinematic lower and upper distance limits, D\ ow and D up 
respectively; and their uncertainties, o"i ow and a up , to conventional radio pulsars as well 
as to otherwise radio-quiet neutron stars with radio bursts; and to supernova remnants 
(SNRs) securely associated with various kinds of neutron stars. The values shown are from 
the stated authors' work, unless otherwise indicated in the table. If the original authors 
gave distance limits meeting our criteria, but neglected to derive unce rtaintie s on these 



FW90. In such 



limits, then we do so ourselves according to the procedure laid out in 
cases, the table entry's reference shows a superscript a. If we judge that the listed authors' 
distance limits themselves need adjustment, we do so and mark the entry's reference with a 
superscript b, and describe details of any such changes in §2.21 Two sources, marked with 
the superscript /, had two or more original sets of distance limits because the cited authors 
evaluated multiple Galactic rotation models without expressing a clear preference for one; in 
those cases we select the one using the standard flat rotation curve, for overall consistency. 
(Note that in none of these cases the various rotation models provided significantly different 
results.) 
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2.2. Notes on individual sources 

In the subsections below, we explain any adjustments to criteria that led to the 
originally published upper and/or lower distance limits. The values themselves are 
summarized in Table [TJ 



2.2.1. SNR Kesteven 73 and AXP 1E1841-045 



Tian &: Leahyl (l2008al ) performed an H I kinematic distance study of the supernova 
remnant Kes 73, which is associated with AXP 1E1841-045. The authors showed that 
the SNR absorption extends to the tangent point, 7.5 kpc distant, which marks the 
lower distance limit. They also made the case that the lack of absorption at v = 84 
km s _1 sets an upper distance limit on the far side of the tangent point. We find these 



arguments compelling, t 



D r 



-up — 

curve. 



10.2 kpc, whereas 



owever, we fin d that t he flat rotation curve then indicates that 



Tian &: Leahyl ( l2008af ) quoted D up = 9.8 kpc for a flat rotation 



2.2.2. PWN G54. 1+0.3 and PSR J1930+1852 



Leahy et al. 



( 120081 ) analyzed H I spectra of PWN G54.1+0.3, which is associated 
with PSR J1930+1852. We confirm that the lower distance limit is at the tangent point. 
While they place the upper distance limit at the Solar Circle on the far side of the Galaxy 
du e to a lack of any negative velocity absorption, we instead adhere to the procedure 



of IFW90I . relaxing the limit to the distance corresponding to the first strong emission 
at negative velocities not showing absorption, i.e., at v — —30 km s _1 . After resetting 
the rotation curve to the flat model with IAU galactic constants, we then find that 
(Aow^up) = (5.0, 12.6) kpc. 
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2.2.3. SNR CTB 80 and PSR B 195 1+32 



Strom Stappera ( 120001 ) measured the H I absorption spectrum of SNR CTB 80, 
which is associated with PSR B1951+32. There is significant absorption out to the tangent 
point, yielding D\ ow = 3.1 kpc. Unfortunately, the published absorption spectrum does not 



extend below v = —15 km s 1 , which is insufficient to establish a _D up measurement. 



Lutz-Kelker Bias and Corrections 



Lutz fc Kelkerl ( 119731 ) first presented the argument that because of the non-linearity of 



sample volume with distance, objects are statistically more likely to be further away rather 



than closer by. Correction for this bias (which is re 
correctly named Lutz-Kelker bias, as discussed by 



ated to the Malmquist bias but is more 



VLMlOl ) is relatively straightforward 



through a likelihood analysis that incorporates probabilities derived from a variety of 



possible measurements. Our derivation is sim ilar to t hat of 
fundamental areas. First, the primary focus of 



VLM10 



but differs in a few 



VLMlOl was biases in parallax measurements, 



while our analysis considers both parallax and distance, which is a more natural quantity 
when dealing with H I distance limits. (Note that the conversion between parallax and 
distance is not a simple inversion in the case of finite uncert ainties, as the transformation 



VLMlOl applied a Bayesian 



between these two quantities is non-linear.) Second, where 
analysis with prior information based on the pulsar luminosity and position in the Galaxy, 
we consider these quantities as measurements and have hence no need for prior information 
at all, which removes the Bayesian character of this analysis and leaves a straightforward 
likelihood analysis. Effectively this is no more than an aesthetic difference, however, which 
does not affect the results. Indeed, our approach could be considered Bayesian with a 
uniform prior. In particular, our analysis considers the following possible measurements: 
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• a parallax measurement, ro meas ; 

• a lower H I distance limit, D\ ow ; 

• an upper H I distance limit, D up ; 

• the pulsar radio flux, S (measured at or near an observing frequency of 1.4 GHz); 

• and the pulsar's Galactic position, G\, G^. 

Given a subset or all of these measurements and assuming no correlations between these 
values, we can determine the probability density function of the pulsar distance, D, through 



p (D\zu meas , Aow, D up , S, G h G h ) = p (D\w meas )p(D\D low )p(D\D up )p(D\S)p(D\G h G h ) . 



In the above equation (as in all equations throughout this paper), we only explicitly state 
dependence on parameters, while dependence on the uncertainties of said parameters is 
implied. In other words, where we write p(D\w meas ), we really mean p(D\zu meas , 0" romeas )- In 
the following, these five terms will be derived; they will respectively be referred to as the 
parallax term, the lower H I limit term, the upper H I limit term, the luminosity term and 
the volumetric or Galactic term. 

3.1. The Parallax Term, p (D\zu mcas ) 

Given a measurement zu me3S with uncertainty cr ro and assuming a Gaussian uncertainty 
distribution, the probability of the true parallax given the data is 



(1) 





Since p (D) = \dw/dD\p (w) oc p (zu) /D 2 , this means 




(3) 
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In the case of asy mmetric uncertainties on parallax measurements (as given, e.g., by 



Chatterj ee et al. 



2009), we assume 



p (za\m meas ) oc H {w 

^meas) exp 



1 / W-n 



(T 



roup 



+H (w meas - w) exp 



1 / ^mcas ^ 



(4) 



with w + o" roup and vj — o" TO i ow respectively the upper and lower limit of the 1 a interval of 
the measurement's probability density function; and with H[x) the Heaviside step function, 
for which 



if x < 0, 
H(x) = l 05 if x = o, 

1 if x > 0. 

For distance, as in the symmetric case, the extra factor of D~ 2 is added, resulting in 



(5) 



p (D\zu mcas ) ^jy^H (1/D - zu mcas ) exp 



1 / Wr, 



1/D 



a 



roup 



+ -^H O meas - exp 



1 / W T , 



1/D 



(6) 



3.2. The H I Distance Limit Terms, p (D\D up , Aowj 



Assuming the distance of the furthest absorbing H I gas is determined to be D\ ow with 
measurement uncertainty cr low , then the probability distribution of the actual distance of 
the limiting gas is given (assuming Gaussian uncertainties) by 

2" 



p(d\D low ) oc 



1 



exp 



low 



1 / .Plow - d 

2 V °"low 



(7) 



where d is the actual distance of the gas and hence the actual lower limit on the pulsar 
distance. This implies that for any pulsar distance D we must have D > d. Hence we derive 



the probability distribution for the pulsar's distance as 

POD 

p(D\D law )= / p(D\d)p(d\D low )dd, 
Jo 

in which 

p(D\d) (xH(D-d) 
with H(x) the heaviside function, as defined above. We therefore have 

POD pD 

p(£>|Aow)oc / H(D-d)p(d\D low )dd = p(d\D low )dd, 
Jo Jo 



(8) 



(9) 



(10) 



which results in 



p(D\D low ) oc 



erf 



/ Aow \ £ I 1 'low 



low 



— erf 



D ]nw - D 



:n) 



with erf(x) = ^ Jq e * 2 dt the error function. 



Analogous to the above derivation, we have the probability distribution for the distance 
of the nearest gas not seen in absorption 

2" 



1 



p(d\D up ) oc — = exp 



27TO". 



up 



1 fD up -d 



2 V <J 



up 



(12) 



which is used in the probability distribution for the pulsar's distance as: 

POD 

p(D\D up )= / p(D\d)p(d\D up )dd 
Jo 



with 



hence: 



p(D\d) xH(d-D), 



p(D\D up )(x / H(d-D)p(d\D up )dd= / p(d\D up )dd 

'0 



which results in 



p(D\D up ) oc 



erf 



V2a 



up 



+ 1 



(13) 



(14) 



(15) 



(16) 
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3.3. The Galactic ("Volumetric") Term, p(D\G h G h ) 



As derived by lLorimer et al.l ( 120061 ). the distribution of pulsars in the Galaxy is not 
homogeneous, but rather follows a distribution of the form 



p(R,ip,z) 



N 
V 



cx R B exp 



If! C ^-R 

E R 



kpc 



(17) 



with N the number of pulsars per volume V and constants Ro = 8.5 kpc, B = 1.9, C = 5 
and E = 330 pc for common pulsars and E = 500 pc for millise cond pulsars (constants from 



model fit C and equations 10 and 11 from 



Lorimer et al. 



Since the volume density is invariant with the coordinate system used, we can use an 
Earth-based coordinate system based on the Galactic coordinates of the pulsar and its 
distance to the Earth, (D, G b , G\), for which p (D, G^, G{) = p (R, ip, z). For the Earth-based 
observer the infinitesimal sample volume now becomes 

SV = D 2 5D5VL (18) 

for a pulsar at given distance D and an infinitesimal solid angle 5fl. The number of pulsars 
in this volume is, hence, 

5N = p (D, Gb, G\) D 2 5D5Vl. (19) 
Since the infinitessimal probability 5P scales with 5N, we get 

p( J D|G b ,G 1 )ocp( J D,G b ,G 1 ) J D 2 . (20) 



If! 5 ^-r 

E R 



Consequently, we derive 

p(D\G h ,Gi) oc R 19 exp 
with 

z(D,G b ) = DsinG h 

and 

R (D, G h , Gi) = Jr 2 + {D cos G h f - 2R D cos G h cos G h 



(21) 



(22) 



(23) 
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3.4. The Pulsar Luminosity Term, p(D\S) 



Finally, since the radio flux, S, of pulsars is related to the luminosit 



of the pulsar 



through SD 2 = L, this measure can be used to constr ain the pulsar distance, t 



i rough 



Faucher-Giguere fe Kaspi 



( 120061 ) 



the luminosity distribution of radio pulsars derived by 
Considering pulsar luminosities at 1.4 GHz observing frequency with luminosity expressed 
in units of mJykpc 2 , they proposed a log-normal function with mean (A) = (log(L)) = —1.1 
and standard deviation cr A = 0.9: 



p (A) oc exp 



With A = log L = log S + 2 log D, we get 



1 / A+l.l 

2 I 0.9 



p{D) oc 



or, given S, 



dD 



1 



p (A) oc — exp 



1 / A + l.l 

2 V 0.9 



p(D\S) oc exp 



1 /logS + 21og J D + 1.1 

2 1 09 



(24) 



(25) 



(26) 



Note that this probability is based on the measured radio flux S of the pulsar, not on the 
H I flux or the luminosity of an associated supernova remnant or the like. Also, given the 
analysis of iFaucher-Giguere fc Kaspil (120061 ) who derived the luminosity distribution that 



we use, the above analysis does not hold for non-radio or bursting pulsars. 



3.5. Combined Distance Probability 



Combining equations [3J CEU [HI EH and into Equation [Q we obtain the complete 
formula for the pulsar distance given the five measurements listed at the start of this 



2 Note that 



Faucher-Giguere &: Kaspil (120061 ) define a "pseudo-luminosity" L = SD 2 that 



avoids the complexities of emission beam and viewing geometries. This approach is practical 
for our purposes and hence we copy their usage of L as an effective 'pseudo- luminosity'. 



section: 



p(D\w, 



Aow, -D up , S, G\, Gb) oc 



1 



1 / ZU- t 



l/D 



) 



2 



meas 





X — 



(27) 



with R as in Equation [221 Note that the parallax term should be replaced by Equation [H] in 
case of asymmetric uncertainties. (Technically Equation [6] can be applied generally to both 
asymmetric and symmetric cases, but for reasons of clarity, we present the more common, 
simplified formula here.) Note also that in case measurements are not available, the relevant 
terms should be omitted, as p (D\m meauS ) (for example) is nonsensical in the absence of a 
zu meas measurement. 
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Equivalently, we find for the pulsar's parallax 

dD ' 



V (^l^meas, Aow, Alp, 5", G\, G h ) OC 



dzu 
oc exp 

1 

x 



P (D\ OTmeas, Aow, -Ap, S, G\, G D ) 
2" 



1 / W r 



zu 



X 



2w 2 
1 

2^ 



2 

erf 



( Aow 



1 +(Tf ' ^ 



ow 

D 



— erf 

- l/w 



D 



low 



1/ 



— 



a/2(Ji c 



y/2, 



a 



up 



X 



1.9 



exp 



IsinGJ _R — Rq 
- 5- 



x — exp 

zu 



wE R 

1 Aog^-21ogcc7 + l.l 

2 1 0^9 



(2f 



where each term contributes a factor w~ 2 , since p(zu\zu meas )5w = p(D\zu mcas )5D, 

p(zu\ D\ ow )Szu = p(D\D\ ow )5D etc., implying that each of the five terms contributes a 5D/Szu 

term. 



VLM10. 



Equation [27] presents the analytic result to the question first discussed by 
However, because in that previous paper parts of the analysis were performed by Monte- 
Carlo si mulation, our present results are more accurate; and in contrast to the analysis by 



VLMlOl . which only considered parallax, we now derive the full formulae for both distance 



and parallax. We therefore present in Table [2] the bias-corrected p arallax a nd distance 



values for the pulsars with parallax measurements first collat ed by 



VLM10 



Results for 



pulsars with H I distance limits (which were not included in IVLMlOl ). are presented in 
Table [TJ For the data in these tables, we use the following definitions: the corrected 
distance (-Dcorr) is the distance for which Equation [27J reaches a maximum; for corrected 
parallax (c^corr) the same convention is used, based on Equ ation 1281 The 1 a uncertainty 



intervals are defined (consistent with IChatterjee et al 



2009) as the narrowest interval that 



contains 68% of the integrated probability density. In practice this means that a level P* 
is found so that the integral of p(D) for those values of D where p(D) > P* , contains 68% 
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of the total probability. For bimodal distributions (which only occur towards the Galactic 
centre, and particularly for the measurements for PSR J1752— 2806), this may result in 
two separate regions (a global optimum and a secondary optimum) which in combination 
contain 68% probability. Estimation of these quantities is analytically unfeasible and is 
therefore performed numerically. The code used to calculate the bias-corrected parallax 
and distance values and uncertainties listed in the tables, is available as a supplement to 



this paper and through an on-line interface on http://psrpop.phys.wvu.edu/LKbias An 
example of the graphical output, showing all five probability terms for PSR J1939+2134 
(B1937+21), is shown in Figure HJ 



4. Discussion and Conclusions 

Of the 80 pulsars with H I distance limits, all but one have post-correction distances 
consistent (at the 1 a level assuming the uncertainties derived from our analysis) with 
the H I limits published. The exception is PSR J2018+2839 (PSR B20 16+28), which has 
a lower distance lim it of 3.2 ± 2.1 kpc, but a parallax measurement of 1.03 ± O.lOmas 



( iBrisken et al. 



20021 ) . which dominates the result and therefore makes the H I distance 
limit irrelevant. Furthermore, there is a single source that is beyond the upper H I distance 
limit (though within la): this is XTE J1810-197, for which we determine a bias-corrected 
distance of 3.7 ± 0.5 kpc, which is just beyond the upper distance limit of 3.4 ± 0.6 kpc 
derived from H I observations. Since for this neutron star both the lower and the upper limit 
are equal; and because no radio luminosity is available, the volumetric term determines the 
slightly higher distance. For 20 sources, the bias-corrected distance is closer than the lower 
H I distance limit (though within 1 a) and 59 (or three out of four) sources are completely 
within the distance limits, with typically bias-corrected distances close to the lower H I 
distance limit. The fact that our analysis finds sources are more likely to be closer to the 
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Parallax (mas) 



Fig. 1. — Example output from our likelihood analysis. For PSR J1939+2134, we show the 
peak-normalised probability distributions of the volumetric (dotted) and luminosity (triple- 
dot-dashed) terms, as well as the distance limits from H I estim ates (dashed and dot-dashed), 



the parallax measurement published by 



Verbiest et al. 



(12003 . thin full line) and the final 



probability distribution for the pulsar distance (thick full line), with peak and 1 a uncertainty 
interval indicated by the vertical lines. The top figure shows these distributions as a function 
of distance, while the bottom figure shows the same distributions as a function of parallax. 
Note that, because of the non-linear relationship between parallax and distance, the most 
likely distance is not necessarily equal to the inverse of the most likely parallax, although 
these values do converge for small uncertainties. 
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lower rather than upper H I distance limit (or, indeed, cl oser even than t he low er limit), is 



Lutz fc Kelkerl ( 119731 ). There are 



unexpected when seen from the perspective presented by 
two reasons for this. 

First, the upper H I distance limits are mostly past the tangent point. This means that 
the volumetric term peaks within - or close to - the range allowed by the H I limits, which 
causes the volumetric bias to be either very weak or non-existent. Second, the pulsars to 
which H I distance limits have been measured, are mostly bright sources, with the exception 
of the flaring neutron stars and those neutron stars that have H I limits derived from 
associations with supernova remnants. The brightness of these pulsars implies a luminosity 
term that peaks at very small distances. 

Comparing the results in Table [T] and the discussion above with the results in Table |2J 
it is clear that the types of neutron star distance estimates (parallax and H I measurements) 
suffer from different statistical biases, although the magnitude of the biases is limited in 
both cases. While parallax measurements are typically biased towards smaller distances 
(i.e. the sources are actually further away than suggested by the measurement) because of 
the relatively limited distance to which this technique works (and the consequential strong 
effect of the volumetric term), the H I measurements are typically biased towards larger 
distances (i.e. the sources are often closer than suggested by the measurement) because the 
volumetric term has little impact and the luminosity term dominates the analysis. 

Finally, of the eight pulsars with both H I distance limits and parallax distances, 
only PSR J1857+0943 (B1855+09) has a bias-corrected parallax t hat is mconsistent w ith 



the parallax measurement. The published value of 1.1 ± 0.2 mas (jVerbiest et al. 



2009|) is 



found to be considerably larger than the most likely value of 0.6ln ? mas, which is partly 



because of the volumetric information (as already found by 



VLM10 



who derived a value of 



0.9 ± 0.2 mas), but also because of the H I limits, which place the pulsar well beyond lkpc. 
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Table 1. Pulsar H I distance limits and Lutz-Kelker-bias corrected distances and 

parallaxes. 



to 

I 



Pulsar name Association ro mcas flow D up Suoo t^Corr -Dcorr Rcf 5 

J2000 B1950 (mas) (kpc) (kpc) (mJy) (mas) (kpc) 

J0141+6009 B0138+59 - 2.6 ±0.7 2.9 ±0.7 4.5 0.30±g;g| 2.3 ± 0.7 (1) 

J0332+5434 B0329+54 0.94 ±0.11 1.7 ±0.7 2.0 ±0.8 203 0.8 ±0.1 1.0 ±0.1 (1,21) 

J0358+5413 B0355+54 0.91 ±0.16 1.4 ±0.7 2.2 ± 0.9 23 0.7 ±0.2 1.0±g;? (1,22) 

J0738-4042 B0736-40 - 2.1 ±0.6 - 80 0.3 ±0.1 1.6 ±0.8 (2) 

J0742-2822 B0740-28 - 2.0 ±0.6 6.9 ±0.8 15 0.16±g;g| 2.0±J;g (3) 

J0837-4135 B0835-41 - 1.8 ±0.8 6.0 ±0.7 16 0.18±g;g| 1-5±J;| (2) 

J0908-4913 B0906-49 - 2.4 ±1.6 6.7 ±0.7 10.0 0.16±g;g| 1.0±q 7 7 (3) 

J0942-5552 B0940-55 - - 7.5 ± 0.7 10.0 0.16^;^ 0.3^ (2) 

J1001-5507 B0959-54 - - 6.9 ±0.7 6.3 0.16±g:g| O.StH (3) 

J1048-5832 B1046-58 - 2.5 ±0.5 5.6 ±0.8 6.5 0.18±g;g| 2.9±l% (2) 

J1056-6258 B1054-62 - 2.5 ±0.5 2.9 ±0.5 21 0.33±g;g| 2.4 ±0.5 (3) 

J1124-5916 SNR G292.0+1.8 core - 3.2 ± 2.0 - 0.08 0.08^;^ 5jJ (4) a 

J1141-6545 - 3.7 ±1.7 - 3.3 0.12±gg| 3±2 (5) a 

J1157-6224 B1154-62 - 3.8 ± 1.4 9.0 ± 0.6 5.9 0.12±ggf 4±2 (2) 

J1224-6407 B1221-63 - 4.3 ± 1.4 11.4 ± 0.7 3.9 0.10±g;gf 4±2 (2) 

J1243-6423 B1240-64 - 4.5 ±1.9 11.5 ±0.7 13 0.10±g;gf 2±2 (1) 

J1326-5859 B1323-58 - 3.0 ±1.0 - 9.9 0.12^°;^ 3±\ (6) 

J1327-6222 B1323-62 - 5.1 ± 1.7 11.8 ±0.6 16.0 0.093tg;^g 4±2 (1) 



Table 1 — Continued 



Pulsar name 
J2000 B1950 



Association 



^mcas 

(mas) 



J1359- 


-6038 


B1356- 


60 


J1401- 


6357 


B1358- 


63 


J1453- 


-6413 


B1449- 


•64 


J1513- 


-5908 


B1509- 


•58 


J1559- 


-4438 


B1556- 


•44 


J1600- 


-5044 


B1557- 


50 


J1602- 


-5100 


B1558- 


50 


J1644- 


-4559 


B1641- 


■45 


J1651- 


4246 


B1648- 


•42 


J1707- 


4053 


B1703- 


•40 


J1709- 


4429 


B1706- 


44 


J1721- 


3532 


B1718- 


35 


J1740- 


-3015 


B1737- 


30 


J1745- 


-3040 


B1742- 


30 


J1752- 


-2806 


B1749- 


-28 


J1801- 


2304 


B1758- 


23 


J1803- 


-2137 


B1800- 


•21 


J1807- 


-0847 


B1804- 


•08 



-Dlow 


D up 


$1400 




-Dcorr 


Ref° 


(kpc) 


(kpc) 


(mjy) 


(mas) 


(kpc) 




5.6 ± 1.7 


_ 


7.6 


09+ 04 

u - uy -0.03 


5±2 


(1) 


1.6 ±0.5 


2.7±0.7 


6.2 


n o-i +0.09 

u - oi -0.06 


1 8+ - 7 


(2) 


2.5 ±0.5 


- 


14.0 


1 3+ 009 


o + l-3 
Z - S -0.8 


(3) 


3.8 ±0.5 


6.6 ± 1.4 


0.94 


u - 1 *-0.02 




(7) 


2.0 ±0.5 


_ 


40 


0-32±g;£ 


9 q+0.5 
z - J -0.3 


(3, 23) 


6.4 ±0.5 


18.2 ± 1.2 


17.0 


O .o8ir 3 2 




(8) 


7.4 ±0.5 


9.4 ±0.4 


5.7 


n n o+0.013 
u - iio -0.007 


8.018:? 


(8) 


4.2 ±0.3 


5.0 ±0.3 


310 


0.21 ±0.02 


4.5 ±0.4 


(1) 


4.8 ±0.3 


- 


16.0 


o.o8ig : g 


5 2+ 2 ' 1 

°- z -0.6 


(9) 


3.8 ±0.5 


_ 


7.2 


0.08lg;g 




(9) 


2.4 ±0.6 


3.2 ±0.4 


7.3 


0-3l±8;8S 


2.6t° :6 


(3) 


4.4 ±0.5 


5.2 ±0.6 


11.0 


0.19 ±0.02 


4.6 ±0.6 


(9) 




5.5 ±0.6 


6.4 


n 20+ - 07 

u - zu -0.03 




(8) 




5.5 ±0.6 


13.0 


n 20+ 008 

U.ZU_g Q 3 


0.2±J"i 


(8) 


0.125 ±0.025 




18.0 


0.08 ± 0.03 d 


2+ 1 - 1 

u - z -0.1 


(1, 10) 


3.5 ±0.9 


6.9 ±0.1 


2.2 


u.iiy_ 005 


4± 1 


(11) 


4.0 ±0.6 


4.9 ±0.3 


7.6 


0.21 ±0.02 


4.4+°;* 


(1) 


1.5 ±0.7 




15.0 


U.ll_ .03 


1 -°-0.9 


(1) 



to 



Table 1 — Continued 



Pulsar name 
J2000 B1950 



Association tn mcas D low D up 5*1400 ^Corr -D Co rr 

(mas) (kpc) (kpc) (mJy) (mas) (kpc) 



Rcf° 



0.5: 

0.13 



0.5 



SGR 1806-20 radiofiarc 2005 - 6.2 ±0.1 - - 0.06 

XTEJ1810-197 radiofiarc 2006 - 3.4 ±0.6 3.4 ±0.6 

J1820-0427 B1818-04 - - 1.6 ±0.5 6.1 

J1823+0550 B1821+05 1.6 ±0.5 - 1.7 

J1824-1945 B1821-19 - 3.2 ± 0.5 - 4.9 

J1825-0935 B1822-09 - - 1.9 ±0.4 12.0 

J1832-0827 B1829-08 - 4.7 ±0.3 5.8 ±0.3 2.1 

J1833-0827 B1830-08 - 4.0 ±0.4 5.3 ±0.3 3.6 

J1833-1034 SNR G21.5-0.9 - 4.0 ±0.3 4.1 ± 0.3 0.071 

AXP 1E1841-045 SNR Kes73 - 7.5 ±1.0 10.2 ±0.3 

J1846-0258 SNRKes75 5.5 ± 0.4 5.9 ±0.5 

J1848-0123 B1845-01 - 4.2 ±0.4 4.8 ±0.4 8.6 

J1852+0031 B1849+00 - 7.1 ± 1.2 16.6 ± 0.9 2.2 

J1857+0212 B1855+02 6.9 ±1.3 - 1.6 

J1857+0943 B1855+09 1.1 ±0.2 1.6 ±0.5 2.0 ±0.4 5 

J1901+0331 B1859+03 6.8 ± 1.4 15.1 ± 0.7 4.2 

J1901+0716 B1859+07 - 2.8 ±0.5 4.7 ±0.8 0.9 

J1902+0556 B1900+05 - 3.1 ±0.4 4.3 ±0.5 1.2 



+0.02 
-0.01 



25+ - 04 

u -^°-0.03 



0.1 

+0.09 
-0.04 



13+4 
3.6 ±0.5 

h0.6 
-0.2 



o.3: 



09+ 07 

u - uy -0.02 



2 0+ 1 - 3 

z - u -0.8 

q 7+I.6 
°- ' -0.9 



+0.2 
-0.1 



0.3 



+0.7 

-0.2 



cist™? 



r 9 +0.5 



0.20 ±0.02 4.5 ±0.5 
0.24 ±0.02 4.1 ±0.3 



102+ 012 
u.iuz_ 005 

0.17±0.01 



9.6±;-« 

r n+0.5 

O-8-0.4 



n7n+ 025 
u.u/u_ 009 



0.21 ±0.02 4.4 ±0.4 
8±2 
8±2 
0.9 ±0.2 
7± 2 



0.08 
0.6 
0.075 



0.20 
0.23 



+0.03 
-0.02 

+0.2 
-0.1 

+0.027 
-0.008 

+0.05 
-0.03 

+0.04 
-0.03 



3.4 
3.6 



+0.9 
-0.7 



+0.6 
-0.5 



(12) 

(13) f 
(1) 
(1) 
(8) 
(8) 

(1) 
(9) 

(14) a,9 

(15) b </ 

(16) a .s 

(1) 
(1) 
(1) 

(1, 24) 
(1) 
(1) 
(1) 



Table 1 — Continued 



Pulsar name Association ro mcas Mow D up 5*1400 ^Corr -Dcorr Ref 6 

J2000 B1950 (mas) (kpc) (kpc) (mJy) (mas) (kpc) 



J1902+0615 


B1900+06 


- 


6.5 ± 1.4 


15.8 ±0.8 


1.1 


071 +0 - 024 

U - U ' ^0.007 


7+3 
' -2 


(i) 


J1903+0135 


B1900+01 


- 


2.8 ±0.4 


4.0 ±0.4 


1.1 


26+ - 04 

u -^°-0.03 


q q+0.6 
o -°-0.5 


(i) 


J1906+0641 


B1904+06 


- 


6.5 ±1.5 


14.0 ±0.5 


1.7 


+0.026 
u - u ' ' -0.006 


7±2 


(i) 


J1909+0254 


B1907+02 


- 


3.8 ±0.5 


- 


0.63 


09+ 05 




(17) 


J1909+1102 


B1907+10 


- 


4.3 ±0.6 


6.0 ± 1.6 


1.9 


14+ - 04 




(1) 


J1915+1009 


B1913+10 


- 


6.0 ± 1.5 


14.5 ±0.8 


1.3 


n 077+0.027 

u - u ' ' -0.009 


7±2 


(1) 


J1916+1312 


B1914+13 


- 


4.0 ±0.7 


5.7 ± 1.7 


1.2 


14+ - 04 


4 -°-0.9 


(1) 


J1917+1353 


B1915+13 


- 


4.8 ± 1.0 


5.7 ± 1.7 


1.9 


0.14 ±0.03 


5± 1 


(1) 


J1921+2153 


B1919+21 






2.8 ± 1.2 


6 




3+ ' 8 

U.O_ 2 


(1) 


J1922+2110 


B1920+21 




4.8 ± 1.8 


16.2 ± 1.0 


1.4 


o.o8i8;g 


4±2 


(17) 


J1926+1648 


B1924+16 




5.2 ± 1.8 


14.9 ±0.8 


1.3 


n n75+ 0028 

U.U / O_o 008 


fi +3 
D -2 


(17) 


J1932+1059 


B1929+10 


2.77 ±0.07 




1.6 ±0.5 


36 


9+ 10 


Sl+ 09 

u - J1 -0.06 


(1, 22) 


J1932+2020 


B1929+20 




4.8 ± 1.8 


14.9 ±0.9 


1.2 


76+ 029 

U.U/D_g q 09 


r.+3 
J -2 


(1) 


J1930+1852 


PWN G54.1+0.3 




5.0 ± 1.8 


12.6 ±0.6 


0.06 


0.085lg;S? 


7 +3 
' -2 


(l 8 )b,g 


J1932+2220 


B1930+22 




10.4 ±0.6 


13.7±0.7 


1.2 


0.08llg;S? 


io.9t : 3 8 


(1) 


J1935+1616 


B1933+16 


n 99 +0.08 
u - zz -0.12 


5.2 ± 1.7 




42 


0.1318:25 


O 7+1.3 

°- ' -0.8 


(1, 25) 


J1939+2134 


B1937+21 


0.13 ±0.07 


4.6 ± 1.9 


14.8 ± 0.9 


10 


o.o8i8:g 




(1, 24) 


J1946+1805 


B1944+17 






1.9 ±0.7 


10 




o.3±8:S 


(1) 
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Pulsar name 
J2000 B1950 



Association tiWas Aow -Dup Shoo wcorr -Dcorr 
(mas) (kpc) (kpc) (mJy) (mas) (kpc) 



Ref 



J1952+3252 B1951+32 SNR CTB80 - 3.1 ± 2.0 - 1.0 

J2004+3137 B2002+31 - 7.0 ± 0.7 12.0 ±0.7 1.8 

J2018+2839 B2016+28 1.03 ±0.10 3.2 ±2.1 - 30 

J2022+2854 B2020+28 0.37 ±0.12 3.1 ± 2.1 - 38 

J2113+4644 B2111+46 - 4.3 ± 0.8 6.5 ± 0.7 19 

J2257+5909 B2255±58 - 3.3 ± 0.7 - 9.2 

AXP 1E2259+586 SNR CTB109 - 4.0 ± 0.8 4.0±0.8 C - 

J2321+6024 B2319±60 - 2.6 ± 0.6 - 12 



3±2 



( 19 )b,h 
(1) 

0.9 ±0.1 0.98lg;JJ (1, 21) 



0.11 
0.092 



+0.06 
-0.03 

+0.023 
-0.009 



0.22; 
0.17; 



-0.07 

-0.03 
-0.02 



13+ 008 



0.21 



0.13 



+0.04 
-0.03 

+0.10 
-0.05 



O 1+0.6 
z - i -0.4 

4± 1 
3± 1 
4.1 ±0.7 



2.7 



+ 1.2 
-0.9 



(1, 21) 
(1) 
(1) 
(20) 

(1) 



CO 



a Distancc limit uncertainty is derived by the current authors. 

b Distancc limit and uncertainty are derived by the current authors. See notes in body of paper. 
c Upper limit based on CO emission from a molecular cloud associated with SNR CTB109. 
d PSR J1752-2806 (B1749-28) has a secondary optimum in its parallax, at 0.15io'.o2 mas - 
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Table 2. Pulsar parallax measurements and Lutz-Kelker-bias corrected distances and 
parallaxes for pulsars without H I limits but with parallax measurements only. 



Pulsar name romcas S1400 ^Corr -Dcorr Ref. 

J2000 B1950 (mas) (mjy) (mas) (kpc) 



J0030+0451 






3.3 ±0.9 


0.6 


1 6+ 10 
u - 3 °-0.07 


n 90+0.10 
1.03 ±0.08 


Lommen et al. (2006, 2000) 


J0034- 


0721 


B0031- 


-07 


„ _ , 1 0.08 
93+ 


11 


Chatteriee et al. (2009) 


J0108- 


-1431 






4.2 ± 1.4 


1.0 


L4 -0.7 

0.37 ± 0.04 


n 01 +0.09 
°- 21 -0.05 

Z ' D -0.2 


Dcllcr et al. (2009b) 


J01394 


-5814 


B01364 


-57 


0.37 ± 0.04 


4.6 


Chatteriee et al. (2009) 


J0437- 


-4715 






6.396 ± 0.054 


112 


6.39 ±0.05 


0.156 ±0.001 


Deller et al. (2008) 


J0452- 


-1759 


B0450- 


-18 


0.64t 06 


5.3 


n -r+0.6 


n ,1+0.2 
0- 4 -o.i 


Chatteriee et al. 


J04544 


-5543 


B04504 


-55 


„ „ . _i_n 04. 


13 


_ _ ,_|_n 04 


„ _ ^_i_n 07 

i-i8±S5:o5 


Chatteriee et al. (2009) 


J05384 


-2817 






72+ - 12 
u ' ' -0.09 

0.80 ±0.35 


1.9 


69+ 011 

4+0-3 
u, ^-0.2 

2.8 ±0.4 


1.3 ±0.2 


Chatteriee et al. (2009); Lewandowski et al. (2004) 


J0613- 


-0200 






1.4 


9+ 04 
u.a_ 2 

n o 9 +0.05 
u -° -0.04 


Verbiest et al. (2009) 


J0630- 


-2834 


B0628- 


-28 


3.0 ±0.4 


23 


Deller et al. (2009b) 


J06334 


-1746 






4.0 ± 1.3 


- 


2+ - 5 


25+ 23 
u - zo -0.08 


Fahertv et al. (2007) 








J06594 


-1414 


B06564 


-1 1 


3.47 ±0.36 


3.7 


3.3 ± 0.4 


0.28 ±0.03 


Brisken et al. (2003) 


J0720- 


-3125 






2.77 ± 0.89 




u - z -o.i 
0.80 ± 0.14 


„ . _|_0 ^ 
u '^-0.1 

_i_o 2 

A -- L -o.i 

q 4 +0.2 

n ^ori + 008 

°- 432 io:oo7 


Kaplan et al. (2007) 


J0737- 


-3039A 






0.87 ± 0.14 


1.6 


Deller et al. (2009a); Bureav et al. (2006") 


J07514 


-1807 






1.6 ±0.8 


3.23 


6 +0 - 6 
2.31 ± 0.04 


Nice et al. (2005) 


J08144 


-7429 


B08094 


-74 


2.31 ± 0.04 


10 


Brisken et al. (2002) 


J0820- 


-1350 


B0818- 


-13 


ri r:i + 0.U3 
U ' O1 -0.04 


7 


51 + u - l, 3 
u - oi -0.04 


1.9 ±0.1 


Chatteriee et al. (2009) 


J08264 


-2637 


B08234 


-26 


2.8 ±0.6 


10 


2.4 ±0.6 


n to + 0.08 
u -° -0.05 


Gwinn et al. ( 1986) 


J0835- 


-4510 


B0833- 


-45 


3.5 ±0.2 


1100 


3.5 ±0.2 


0.28 ±0.02 


Dodson et al. (2003); Backer & Fisher (1974) 


J09224 


-0638 


B09194 


-06 


0.82 ±0.13 


4.2 


0.82±g : « 


J - L -o.i 


Chatteriee et al. (2001) 


J09534 


-0755 


B09504 


-08 


3.82 ±0.07 


84 


3.82 ±0.07 


0.261 ± 0.005 


Brisken et al. (2002) 


J10124 


-5307 






1.22 ±0.26 


3 


1.11 ± 0.25 


7 +0 ' 2 

u -'-o.i 


Lazaridis et al. (2009) 


J10224 


-1001 






1.8 ±0.3 


3 


1.7 ± 0.3 


52+ ' 09 
u - oz -0.07 


Verbiest et al. (2009) 


J 1024- 


-0719 






1.9 ±0.4 


0.66 


1.5 ±0.4 


49+ 12 
u, ^ a -0.08 

23+ ' 17 
u - zo -0.07 

0.35 ±0.02 


Hotan et al. (2006) 


J1045- 


-4509 






3.3 ± 1.9 


3 


3 +0 - 4 
2.80 ±0.16 


Verbiest et al. (2009) 


J11364 


-1551 


B11334 


-16 


2.80 ±0.16 


32 


Brisken et al. (2002) 


J12394 


-2453 


B12374 


-25 


1.16 ±0.08 


10 


1.16 ±0.08 


0.84 + 0.06- 0.05 


Brisken et al. (2002) 


J13004 


-1240 


B12574 


-12 


1.3 ±0.4 


2 


1 o+ - 4 


6+ 02 


Wolszczan et al. (2000) 


J1456- 


-6843 


B1451- 


-68 


2.2 ±0.3 


80 


2.1 ± 0.3 


43+O.O6 


Bailes et al. (1990); Manchester et al. (1980) 


J15094 


-5531 


B15084 


-55 


0.47 ± 0.03 


8 


0.47 ± 0.03 


2.1 ± 0.1 


Chatteriee et al. (2009) 


J15374 


-1155 


B15344 


-12 


0.98 ±0.05 


0.6 


0.97 ± 0.05 


1.01 ± 0.05 


Stairs et al. (2002) 
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Table 2 — Continued 



Pulsar name 
J2000 B1950 



tX>meas £1400 t^Corr Dqoii 

(mas) (mJy) (mas) (kpc) 



Rcf. 



J1543+0929 B1541+09 0.13 ±0.02 5.9 

J1600-3053 0.20 ±0.15 3.2 

J1643-1224 2.2 ±0.4 4.8 

J1713+0747 0.94 ±0.05 8 

J1744-1134 2.4 ±0.1 3 

J1856-3754 6.2 ±0.6 

J1900-2600 B1857-26 0.5 ±0.6 13 

J1909-3744 0.79 ±0.02 3 

J2022+5154 B2021+51 0.50 ± 0.07 27 

J2048-1616 B2045-16 1.05^°;^ 13 

J2055+3630 B2053±36 0.17 ±0.03 2.6 

J2124-3358 3.1 ±0.6 1.6 

J2129-5721 1.9 ±0.9 1.4 

J2144-3933 6.05 ± 0.56 0.8 

J2145-0750 1.6 ±0.3 8 

J2157+4017 B2154+40 0.28 ± 0.06 17 

J2313+4253 B2310±42 0.93lg;°f 15 



0.16 ±0.02 5.9] 



021 +o.io 

1.9 ±0.4 
0.93 ±0.05 
2.4 ± 0.1 
6.0 ±0.6 



2.1 



+0.9 



0.42 



+0.09 

-0.06 



-+0.06 
-0.05 



1.051 
0.42 ± 0.02 
0.16" 



,+0.02 
-0.01 



0.3 



+0.3 
-0.1 



0.7 



+0.4 

-0.2 



0.79 ± 0.02 


1.26 ±0.03 " 


0.49 ± 0.07 


1 H+0.3 


1 05+ ' 03 
i ' uo -0.02 


95+ 02 
U - MO -0.03 


0.17 ± 0.03 


5 +0 - 8 
°' u -0.6 


z -°-0.7 




5+ ' 6 
U -°-0.3 


0.41°;? 


5.77 ± 0.57 


lfi+ 02 
u - iD -0.01 


1.5 ±0.3 


57+ - 11 
u - d '-0.08 


0.29 ±0.05 


2 q+0.5 


QO+ ' 06 
u - MZ -0.07 


1 06+ 00S 



Chattcricc ct al. (2009) 
Verbiest et al. (2009) ; Jacobv et al. (2007) 
Verbiest et al. (2009) 
Verbiest et al. (2009) 
Verbiest et al. (2009) 
van Kerkwiik fc Kaplan (2007) 
Fomalont et al. (1999) 
Verbiest et al. (2009) ; Jacobv et al. (2003) 
Brisken et al. (2002) 
Chatteriee et al. (2009) 
Chatteriee et al. (2009) 
Verbiest et al. (2009) 
Verbiest et al. (2009) 
Deller et al. (2009b) 
Verbiest et al. (2009) 
Chatteriee et al. (2009) 
Chatteriee et al. (2009) 



